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Abstract. We consider classical particles on the line with the Weierstrass p 
function as potential. This system parameterizes special solutions of the KP 
equation. We derive the trace formula which relates the Hamiltonian of the par- 
ticle system to the residues of some Abelian differential (meromorphic one-form) 
on the spectral curve. Such formula is important for the construction action-angle 
variables and study invariant Gibbs' states. 



1. Introduction. The subject of this note is a system of classical particles 
on the line interacting with the Hamiltonian 

N 2 ^ 
Ti=l n,m=l 

The parameter* a = 1 corresponds to attractive particles and a = i = 
corresponds to repulsive particles. The potential is the Weierstrass p function 
with real period 2uj and pure imaginary period 2a;'. The system includes the 
well known integrable potentials sin~^ x, sinh~^ x and x"^; which correspond to 
various degeneration of p. 

In a remarkable article Airault, McKean, Moser, [AMM], discovered a con- 
nection between particles with rational or elliptic potential and the Korteweg-de 
Vries equation 
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Few years later Krichever, [K], found an isomorphism between particles with el- 
liptic potentials and special solutions of the Kadomtzev-Petviashvilli (KP) equa- 
tion 

^2 / *^ ^ \ 

— (7 Uyy — {^Uf -\- —UUx '^Uxxxjx- 

The case a = 1 corresponds to the KP-2 equation and cr = z to KP-1. The KP 
equation arises as a compatibility condition for the zero curvature representation 

[ady -L2,dt- Ls] = 0, 

where 

L2 = dl- u, 



X 

-^3 = 9| - ^udx - W. 



Any such solution is associated with a spectral curve Fjv of genus N, defined by 

N 



RN{k, z)^Yl ^n{z)k'' = det(L + 2k), 



n=0 



where L is a N x N matrix which depends on q,p and z. The functions rn{z) 
are elliptic, so the curve Fjv is an A^-sheeted covering of the elliptic torus. The 
matrix L has a simple pole above z = and can be expanded in powers of z 



L= -1^-^^ + 0(1), 



where = — 2(1 — 5^^) is a constant matrix. This "zero order" approximation 
provides all the information needed to solve the direct spectral problem and 
obtain a formula for the solution in terms of Riemann theta functions, see [K]. 
In this note we address a different question. Is there a formula of the type 



TV 



where are parameters of the Riemann surface associated with the system? 
We give an affirmative answer to this question here. The formula is needed, see 
[MCVl-2] , to express the canonical measure as 

e-^dvol = e-" \{dl d(j) = e-^ ^ W dl dcf), 
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where d vol is produced from the basic symplectic structure* u = dp A dq and 
Ps and (/»'s are classical action-angle variables constructed from uj. After this is 
done one can try to compute the partition function, [VI]. 

Now analysis of the direct spectral problem requires a "second order" approx- 
imation 

Z = li:(-i) + L(«) + ^i:W + o(z2), 

z 

Such an approximation provides the coefficients k^^ and k^^ (Theorem 9) for 
the expansion of the function k 

k,{z) = ^ + k^^^ + zk['^ + 0{z') 
z 

on the "upper" sheet of the curve. The desired formula can be easily obtained by 
a simple application of Cauchy's theorem. Moreover, in the repulsive case the 
parameters I'^ are real for all configurations of particles. In order to prove this 
we show that in the expansion 



kaiz) = -- + k^^^ + Oiz), 



a = 2, • • • , AT. 



the coefficients ka^ are distinct for all a = 2, • • • , AT and a generic configuration 
of particles. That much information can be obtained by perturbation techniques 
on "lower" sheets. 

Presumably are moduli of the corresponding A-sheeted covers and the 
actions relative to some symplectic structure u' on the phase space, but this is 
not proved. Compare [V2] for the case of the cubic Schrodinger curves. We will 
return to this issue elsewhere. 

2. Elliptic solutions of KP hierarchy. Consider the particle Hamiltonian 
on the line 

N 2 ^ 
n=l n,m=l 

The Hamiltonian produces a system of first order equations of motion 
9H 

qn= -K—, n = l,...,N, 

OPn 

Pn = = 4cr^ p'{qn -Qm), n=l,...,N. 



*Here uj is not the same as 2w, the period of the elliptic functions. 
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The system can be written in the form 

(1) *g n = 4cr^ ^ p'{qn - qm), n = l,...,N. 

The key step in the embedding of the particle system into the eUiptic solutions 
of the KP equation is the following theorem 

Theorem 1. [K]. The equations 

N 

ady-dl + 2j2p{x-qn{y)) V'^O, 



n=l 

n 



ady -dl + 2Y^ p{x - qn{y)) 



n=l 







have solutions of the form 



N 



ipix, y, k,z)^^ an{y, k, z)^{x - qn, 2;)e'=^+ 



N 



^\x, y, k,z) = Y^ al{y, k, z)^{-x + qn, z)e~ 



■kx—<7 k y 



n=l 



where* ^{x,z) = ^^^^e^*'^''^? if and only if qn{y) satisfy the system of equa- 
tions (1). 

The proof is obtained by requiring that singularities of the form [x — qn)~'^ 
and {x — qn)~^ vanish. This condition can be written in a compact form with 
the aid of N x N matrices L and M 

\ s^n I 

Lemma 2. [K]. The vectors a{y,k,z) and a^{y,k,z) satisfy the equations 

{L + 2k)a = {ady + M)a = 

and 

a\L + 2k) = a\ady + M) = 0. 

These equations determine the curve Tn which is the subject of the next 
section. 



*a{z) denotes the Weierstrass function. 
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3. Riemann surface. The matrix L can be simplified using a gauge transfor- 
mation 

L = GLG-^, 

where G„m = e'^'^^^'^^Snm- Then 



Lnm = CrPnSnm + 2$o(gn " 9m, z){l - 6nm) 

and 



a{z)a{x) 

The existence of a nontrivial vector a : {L + 2k)a = imphes that -Rat (A;, z) = 
det(L + 2k) vanishes and this condition determines the curve. We denote by 
P, Q, etc., points {k, z) on the curve. 

Lemma 3. [K]. The determinant R^ik^z) can he written in the form 

N 

RN{k,z) = Y,rn{z)k'', 

n=0 

where r^iiz) are elliptic functions of z. 

The curve F^v is an sheeted covering of the elliptic curve. The next lemma 
describes symmetries of the curve. 

Lemma 4. (i). a = 1. The curve Fjv admits the antiholomorphic involution 

Ti : {k,z) ik,z). 
(ii). a = i. The curve Fjv admits the antiholomorphic involution 

n- {k,z) {-k,-z). 



Proof, ii). For rectangular lattice a{z) — a{z) and (^q{x^z) — (^q{x^z). There- 
fore RN{k, z) — RN{k, z). 



{ii). Note first ^o{x, -z) = -^q{-x,z) and ip+2k= -{ip + 2{-k)). There- 
fore RN{k,z) = {—l)^RN{—k, —z). The proof is finished. □ 

The function ^q(x,z) has a simple pole at 2; = and can be expanded in 
powers of z 

$o(x, z) = -- + Cix) + lipix) - eix))z + Oiz'). 
z z 
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Therefore*, 

L{p, q, z) = -L(-i) {p, q) + L^") (p, q) + zL^^^ {p, q) + 0{z^), 
z 

where 

L^nm = (^PnSnm + Ki^n " 9m)(l " ^nm), 
-^nm = (P - C^)(9n - 9m)(l - ^nm)- 

The matrix L^~^^ is a constant matrix. Its spectrum and eigenvectors can be 
easily computed. 
Let 



27T 

and /5q, = — (a-1), Q;=l,...,iV. 



Then 

(l^-i) + 2A;a) aa = 

with ka = —1 ioT a = 2, . . . , N and ki — N — 1. 

The necessary information about the curve is obtained by perturbation of this 
trivial case. The eigenvectors aoi{z) and eigenvalues ka{z) can be expanded in 
power series in z 

aa(z) = a^°) + a'^^h + a^^^z^ + . . . , a^°) = a^, 

Z 

The index "ct" labels the sheets of the curve Fn. We call the sheet a = 1 the 
"upper" sheet. In fact, the upper sheet is distinguished by k^~^\ the leading 
term of the asymptotics. The "lower" sheets {a > 2) are distinguished for generic 
configuration of particles by different values of k^^"* since all k^~^^ = —1. This is 
proved in Lemma 5. The proof of Lemma 6 shows that all leading terms are 
distinct. This implies that all "lower" sheets can be indexed according to these 
asymptotics and aa^ = aa- 



*We omit ~ above L and a to simplify the notations. 
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Lemma 5. For generic configuration of particles k^^^ are distinct 

Proof. To prove the statement we need first order perturbation theory for mul- 
tiple eigenvalues. 

We choose — 1 vectors e^°) in the subspace generated by ai°^ , a = 2, . . . , A/"; 

N 
a=2 

where the 77's depend on 6^°^ and are such that 

(^L{z) + 2k{z)^ e^{z) = 0. 
L{z), kj{z), e^{z) can be expanded in integer powers of z 

kfi^z^ = ki^ ^ -\- ky ^ -\- zk^ ^ -|- . . . , 

z 

Now we collect terms in the identity 

_z z 
Terms in z^ produce 

i:(-i)e(i) + L(o)e(o) + 2k^-^^e^^^ + 2A;(°)e(o) = 0, 

and 

(L(-i)e(i),ai°)) + (LWeW,ai°)) + 2/fc(-i)(e(i),aio)) + 2/c(°)(e(°),ai°)) = 0. 
Using the selfadjointness of L^~^^ 

(L(-^)eW,a(°)) = (eW,L(-^)a(°)) = -2^(-^)(eW, a^). 

Therefore 

(LWe(°),aio)) = -2jt(°)(e(«),aio)) 

or 

f;r;«(L(«)ai°),ai«?) = -2/c(«)r;«. 

a=2 

The eigenvalues —2k''^^ of the matrix (L^'^^ai^'', a^^)), ol^o! = 2, . . . ,N; are dis- 
tinct for generic configuration of particles. □ 



= (0) 



= 0. 
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Lemma 6. For all configurations of particles the "zero" order approximations 



of the eigenvectors a^iz) = a^^^ + a^^^ z -\- . . . of the spectral problem {L + 2k)a = 
normalized by the condition a^'^\l) — 1 are given by the formula 



( \ \ 



and l3o, 



27r 
iV 



(of-l), Q;=l,...,iV. 



ygi/3c«(iV-i) J 



Proof. If ^ as n X n matrix then AA^ = det A /, where A^ is the matrix which 
consists of auxihary minors of see also [KNS]. For any column r = 1, . . . , N 

r~ 1 ^ 
L + 2A; 



a{p) 



pr 



L + 2k 



Ir 



in the vicinity of z — 0*. 

We know that k{z) = - - + k^^'> + ... and 



L = --[E-i] + m + ..., 



Then 



Since rank £■ = 1 we have 





A 


L + 2k 






pr 



-— + [L(«) + 2fc(«)] + 
z 



pr 



ipr Z — ; 




where the subscript s means that the s-th column is replaced by 



For a generic configuration k^^^ are distinct on all sheets of the curve. This and 
the formula for a(p) imply that all a^^^ are also distinct and therefore match ai°\ 
Since all ai°^ are fixed the statement is true for all configurations of particles. □ 

The following two lemmas are simple consequences of the discussion above 



*Note that a is not a function on the curve, since the entries of the matrix L + 2fc are not 
elliptic functions. 
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Lemma 7. [K]. The determinant RN{k, z) can be written in the form 

RN{k, z) =2^{k- + +...)) n - (-;+ +■ • ■)) • 

The genus of the curve V n can be easily computed using the Riemann-Hurwitz 
formula. 

Lemma 8. [K, KBBT]. (i) The elliptic case: 2u!,2u' < oo. For generic config- 
uration of particles the genus of the curve F^v is N. 

(ii) The rational case: 2a;, 2a;' = oo. The genus of the curve Fjv is 0. 

4. Asymptotics for k{z). The main result of this section is the following 

Theorem 9. On the "upper" sheet for ki{z) the following asymptotics hold 

k,U)-^-^ + z - ^] + O(z^) 



where 



and 



N 2 ^ 
n=l n,m=l 



N 



n=l 



To prove the theorem we need the following lemma which is the second order 
perturbation theory of simple eigenvalues adapted to our considerations. 

Lemma 10. The following identities hold 



2/c(°) = -(i:(o)a(o),a(°)), 



a=2 

where 

N 



{f,9) = ^^fngn, 



N 

n=l 
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for any f,g eC^. 

Proof. We start with the identity 

z 



+ 



+ 2/c(°) + 2A;(i) + . . .)(a^°) + za^'^ + ...) = 0. 



Collecting the terms in 2; ■'^ 



Note, 



,(0) 



,(0) 



1 



Therefore = iV — 1, as it has to be. 

Now we collect terms with z^ 



(2) 
and 



L(-i)a(i) + L^aW + 2k^-'^a^'^ + 2fcWaW = 



(3) (L(-i)a(i),a(°)) + (L(«)a(«),aW) + {2k^-'^a^'\a'^'^) + (2A:(«)a(°), a(«)) = 0. 
Note 

(4) (L(-i)a(i),a(°)) = (a^^), L(-i)a(o)) = {a^'\-2k^-'^a^'^), 

which implies that the two terms in (3) cancel each other. Prom the definition 
{a^'^\ a^^^) = 1 we obtain the first statement of the lemma. 

Now we derive the formulas for a*-^-* on the "upper" sheet, which will be useful 
later on. Multiply (2) by ai°\ a = 2, . . . , N, 

(5) (L(-^)aW,af ) + (L(°)a(°),ai°)) + (2/c(-^)a(^), ) + {2k^'^ a^'\ a^^^) = 0. 

The last term vanishes due to the orthogonality of eigenvectors corresponding 
different eigenvalues. Similar to (4) 

(i:(-i)a(i),aW) = (a(i),L(-i)aW) = 2(aW,aW) 
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QL(-i)+2(-i))ai°)=0, a = 2,...,N. 

Therefore 

2(a«,ai°)) + (LWa(°),ai°)) + 2(iV- l)(a«,ai°)) = 

and 

{a^'\ai^^) = -^{Li'^a^'\ai^^ a = 2,...,N. 

The condition ||a(2;)|p = 1 + 0(2;^) imphes 

1 = {a{z), a{z)) = (a(°) + za^''^ + ..., a^^) + za^^^ + . . .) 

= (a(°), aW) + ^ [(a(°), a^^)) + (a^^), a(°))] + . . . . 

and (a(^),a(°)) = 0. Therefore 

JV ^ TV 

a=2 a=2 

In order to prove the second formula of the lemma we collect terms with z^ 

L(-i)a(2) + L(o)^ii) ^ ^(1)^(0) ^ 2yfc(-i)a(2) + 2k^'^a^'^ + 2k^'^a^^^ = 

and 

(L(-i)a(2),aW) + (LWa(i),a(o)) + (L^^^a^, + 
(6) + 2/t(-^) (a(2) , a(°)) + 2/c(°) (a^^) , a^^) ) + 2k('^ (a^") , a^")) = 0. 

Note 

(L(-i)a(2),a(0)) = (a(2),L(-i)a(°)) = -2k^-'\a^'\a^'^) 

and two of the terms in the formula (6) vanish. Using the normalization condi- 
tions 

(a(°),a(°)) = 1 and (a(^),a(°)) = we obtain 

2it« = -(L(0)a«,a(°)) - (L«a(°), a^). 
Using the formula for a^^\ we finally have 

a=2 
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The lemma is proved. □ 

Now we use the lemma to compute the coefficients k^^^ and k^^\ Note, first, 
that L(°) can be split into two parts L^^^ = a A + B, where 

The matrix A is symmetric and B is skew-symmetric 

N ■ 

The term with B vanishes due to skew-symmetry. 

The computation of the next term k^^"^ is much more involved. Using the 
decomposition for L^^^ we have 

N 



a=2 

-(LWa(o),a(°)) = / + //. 



Step 1. Our goal is to evaluate the first term /. Arguments for the attractive 
and the repulsive case are different and we consider these two cases separately. 
Repulsive case, a = i. Using the symmetry of A and skew-symmetry of 5, 



N 



-2^f:|z(^a(°),a(°)) + (W°),a(°))|^ 



a=2 

TV 



2N 



(0)\|2 



a=l 



The last term can be easily estimated 

1 



2N 



|(yla(°),a(°))|2 



1 



p2 



To estimate the first term, we introduce a' such that 

a' = N - a + 2. 



TRACE FORMULA FOR THE SYSTEM OF PARTICLES 



13 



Then 



and 



Now 



2ir 2tt 
/3a' = ^(a'-l) = -(iV-« + l) 



/?a + /3a' = (mod27r). 



2N 



i-f;K(Aa(°),a(°)) + (W°),a(«))p 



N 



AN 



a=l 
N 



AN 



Using the identity \a + + \a - 6p = 2|ap + 2|&p we obtain 



Atractive case, a = 1. Again using properties of A and B we have 



a=2 
AT 



a=2 



1 N 
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The second sum vanishes. Indeed, 

N 

2$5 5](^a(0),aL°))(W0),aW) 

AT 



(Aa(0),ai°))(W°),a^)) + {Aa(o),a^^)){Ba^'\a^^)) 



N 



=0. 

Therefore, 

AT 



AT- 



Combining the results for two cases cr^ = ±1, we obtain 

N 

Step 2. Now we estimate £ | (Aa^^), ai°^)|2 and |(5a(°), ai°^)|2. To do this 
we introduce for two polynomials 

AT AT 

fc=i fc=i 

the inner product 

AT 

< P, Q >= ^ P(e^^« ) 

a=l 

It is easy to see that 

AT 
k=l 

Therefore 

AT N ^ N 

5:i(Aa(°),a(°))r = E 1^ E^'n^"^'"^""^!' 

a=l a=l n=l 

^ AT AT ^ N 

n=l n=l n=l 
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Similarly, 

N ^ N N 



Finally, 



AT 

a=l n=l s=l 



2 AT 1 ^ 2 

n=l n=l 

2 1 ■'V 2 

- V ^ p2 



n=l 

^^ep 5. Now using the expression for //* 



we will estimate 



N 

s^s' 



1 ^^^o.. ^V 
n=l \s=l / 

We will prove that 

^ iV / AT \^ 1 

(7) = ^^^]y2^ E ^(^^^ ~^«')- 

This would imply 

(7^ 1 ^2 „ 

and complete the proof of the theorem. 



2 



*In the rational case p — (^^ = 0; L^^^ and // vanish identically. 
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Step 4- To prove the identity (7) we make some manipulations with the sums 



TV 



7V2 



n=l 



N 



CiQn - Qs) 



,s=l 



s^s' 



iV2 



s^s' 



1 

1 



{N-2)J2CHQs-qs')-'2 Yl Ci'in-qs)C{qn-qs 



Y -qs') + 2C(9s - qn)C{qn - qs') 



Fix ki < k2 < and consider all their permutations. Collect the terms with 
such indices 



2 



Y (^'^i - 9fc2 ) + {qk2 - 9fe3 ) + {qk3 - qki ) 

ki<k2<k3 

+2C(9fci - 9fc3)C(?fc3 - ?fc2) + 2C(?fci - qk2)C{qk2 - qka) 

+2C(9fc2 - qki)C{qki - qka) 

2 

= 7^2 (C(9fci -9^2) + Cfes -9/03) + C(%3 -9fei))^ 

ki<k2<k3 

Using the identity 

(C(«) + C(^^) + C{s)f = p{u) + p{v) + p{s) 
for u + V + s = 0, we obtain 

jp ^(^fei -ik2) + p{qk2 - ?fe3 ) + p{qk3 -qki)- 

ki <k2<k3 

The last expression is invariant under permutations and therefore it is equal 
2 

Y piis - 1^') + piis - qn) + piqs' - qn) 

N-2 



N2 



Y -(Is')- 



s^s' 
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This completes the proof of the identity (7) and the theorem. 

Example. For N = 2 the equation defining the curve T2 has the form 



2 . ,(7P /a2p2 ^2jj 



k' + k—+{— —-p{z)]=0. 



Solving the quadratic equation 



Expanding k\{^z^ at 2; = we obtain 

5. Trace formula. If the total momentum vanishes, Pjv = 0, then the result 
of Theorem 9 becomes 

(8) h{z) = + z^^ + O^z"). 

Theorem 11. The following identity holds 

2 ^ 

where I'^ = —ka^. 

a = i. The antiinvolution Ti does not permute the sheets of the curve. The 
variables I'^ are real for all configurations of particles. 

a = 1. For some configurations of particles the antiinvolution ti permutes 
some lower sheets a and r\a. It also leaves the other sheets invariant. Cor- 
responding variables I'^ and 1'^^^ form complex conjugate paires. All other I's 
corresponding invariant sheets are real. 

Proof. By Cauchy's theorem 

TV 

k{P)p{z{P))dz{P) = -J2^J^ k{P)p{z{P))dz{P), 

where 70 is a small contour surrounding Pa, the point on the cc'th sheet above 
z = 0. The asymptotics (8) implies the result. 
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The second part of the proof is different for repulsive and attractive cases. 
We treat them separately. 

Repulsive case, a = i. If the pair [k, z) satisfies 

/,(,) = i4-i) + /,(0) + 4i), + ... 
z 

in the vicinity of Pq., then the pair [—k,—z) also satisfies similar expression in 
the vicinity of Prio: 

— k{z) ~ ^Z^^TiOi^ + ^r°a + ^rla(~^) + 

Therefore 

z 

Comparision shows that /ci"'' = }St% odd and k^a^ = —ki^l for n even. 

We know that k[~^^ = N — 1 and = — 1 for ck = 2, . . . , A?" so that Tj 

leaves the upper sheet invariant. 

The Lemma 5 states that the values of ka^ are distinct on different sheets 
of the curve for a generic configuration of perticles. The matrix {L^^^aa\a^^)) 
(see the proof of Lemma 5) is skew symmetric and therefore all ka^ are pure 
imaginary: ka^ = —ka^. These together with the identity ka^ = —k^X imply 
that TjO; = a and antiinvolution Tj does not permute sheets. Therefore ka^ = ka^ 
and approximation arguments complete the proof. 

Attractive case, a = 1. Arguments as before lead to k^^ = k^a for all n. 
Skew-symmetry of the matrix L^^^ is lost but for generic configuration of particles 
ka^ are distinct. The antiinvolution ri leaves the upper sheet invariant, but is 
can permute the lower sheets of the curve. An example after the proof of the 
Theorem shows how it happens for N = 3. Therefore, some lower sheets a and 
Tia are interchanged by antiinvolution ti and some are invariant. These and 
fci = ^Tia\ together with approximation arguments complete the proof of 
Theorem. □ 

Example. For N = 3we can explicitly compute expansion for ka{z) at ^ = 0. 
In this case 

Rsik, z) = 8k^ + AaPk'^ + [a'^{P'^ - 2H) - 2Ap{z)\ k + 8aJ + 8p'(z), 
where 



8 J = a'^pip2P3 + Apip{q2 -qz)-\- 4p2p(?i - 93) + 4p3p(9i - ?2)- 
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If P = the equation of the curve becomes 



k-'-k 



+ aJ+p'{z) = 0. 



Cardano's formula for the roots 



k^+pk-q = 



has the form 



ki{z) 



ks (z) = u + + ^|-v^, 



where uj 
obtain 



e 3 



R 



+ 



27' 



After elementary but tiresome calculations we 



ki{z)^ -+z^+0{z') 
z lo 



k2{z) = — + —— + z 
z 6 



J 



k3{z) = 



1 aVSH 



6 



+ z 



J 



+ Oiz^), 



+ 0{z'). 



Consider the case (7=1. If H < 0, then ti permutes lower sheets of the 
curve Fs; I2 and Is are complex conjugate. If H > 0, then ti leaves lower sheets 
invariant; I2 and Is are real. 

Lemma 12. Let Pn = 0, n = 1, . . . , N. Then the curve Fjv admits involution 

T_ : (k, z) {-k, -z) 

The involution r_ leaves the upper sheet invariant and k[^'^ = for N even. It 
can permute lower sheets and in this case k^^ — k^a for n odd and k^^ = 
—k^a for n even. 



1. All k 



(0) 



a 



N are pure imaginary. 



Proof. Using the identity $(g, —z) — — $(— g, z) we have 



L{q,p = 0, -z) + 2{-k) = [L{q,p= 0, z) + (-1)^. 
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Therefore R{—k, —z) = {—l)^R{k, z) and the existence of t_ is proved. 
If a pair (A;, z) satisfies 

/,= l^(-i) + ;,(o) + ^(i),+ .... 
z 

Then ^ 

— k = -kr_a + k^]a + ^^-a{~^) + 

—Z 

Therefore k^^ — k^Ja for n odd and /ci"''* = —k^Ja for n even. Since k\ = 

N — 1 and ka^ = — 1, ck = 2, . . . , A'"; t_ leaves the upper sheet invariant. The 
involution r_ can permute the lower sheets. It is demonstrated by the example 
for N = 3. 

(7 = 1. Skew-selfadjointness of the matrix L*^°) is restored under the conditions 
of the Lemma. □ 

6. Appendix. The Weierstrass a{z) has periods 2a; and 2a;' and is defined as 

where = 2un + 2u'n' and Yl is taken with n, e Z^; + n'^ > 0. ({z) and 
p{z) are difined similarly: 

d d 

C(^) = log ^(^)' p(^) = -:^-C(^)- 

For more information see [HC]. 

Acknowledgments. Finally, I would like to thank I.M. Krichever and H.P. 
McKean for numerous stimulating discussions. 
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